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ABSTRACT 


A buckling analysis was performed on a hat-stiffened panel subjected to shear loading. Both local 
buckling and global buckling were analyzed. The global shear buckling load was found to be several times 
higher than the local shear buckling load. The classical shear buckling theory for a flat plate was found to 
be useful in predicting the local shear buckling load of the hat-stiffened panel, and the predicted local shear 
buckling loads thus obtained compare favorably with the results of finite element analysis. 
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NOMENCLATURE 

2 

cross-sectional area of one corrugation leg, A = lt c , in 

cross-sectional area of global panel segment bounded by p, 

A = A+pt s + l(f l -f 2 )t c , in 2 

Fourier coefficient of assumed trial function for w (x, y), in. 
length of global panel, in. 

horizontal distance between centers of corrugation and curved region, 

b = \ f 2 ) ’ in- 

width of rectangular flat plate segment, in. 
width of global panel, in. 

E s t] 

flexural rigidity of flat plate, D = j-, in-lb 

12 ( i _ v ) 

S 

flexural stiffness parameter, J D x D y , in-lb 

transverse shear stiffnesses in xz-, yz -planes, lb/in 

effective bending stiffnesses of equivalent hat-stiffened panel, in-lb 

one-half of diagonal region of corrugation leg, in. 

2 

modulus of elasticity of hat material, lb/in 

2 

modulus of elasticity of face sheet material, lb/in 
lower flat region of hat stiffener, in. 
upper flat region of hat stiffener, in. 

2 

shear modulus of hat material, lb/in 

2 

shear modulus of face sheet material, lb/in 

distance between middle surfaces of hat top flat region and face sheet, 
h = h c+\(‘c + t s)’ in- 

distance between middle surfaces of hat upper and lower flat regions, in. 



0 corrugation angle (angle between the face sheet and the straight diagonal 

segment of corrugation leg), rad 

v Poisson ratio of hat material 

C 

v Poisson ratio of face sheet material 

S 

2 

x shear stress, lb/in 

xy 

( ) critical value at buckling 


INTRODUCTION 

Recently, various hot-structural panel concepts were advanced for applications to hypersonic aircraft 
structural panels. Among those panels investigated, the hat-stiffened panel (fig. 1) was found to be an ex- 
cellent candidate for potential application to hypersonic aircraft fuselage panels. This type of panel is 
equivalent to a corrugated core sandwich panel with one face sheet removed. 

Buckling behavior of the hat-stiffened panel under compressive loading in the hat-axial direction, was 
investigated by Ko and Jackson recently (ref. 1). They calculated both the local and global (general panel 
instability) compressive buckling loads for the panel. The calculated local compressive buckling load was 
found to be far lower than the global compressive buckling load, and compared fairly well with the exper- 
imental data. To fully understand buckling characteristics of the hat-stiffened panel, the shear buckling be- 
havior of this panel needs to be investigated. 

This report presents the local and global buckling analyses of the hat-stiffened panel subjected to shear 
loading. The predicted shear buckling loads are compared with the finite element shear buckling solutions. 

SHEAR BUCKLING ANALYSIS 

To analyze the buckling behavior of the complex structure shown in figure 1, two approaches were 
taken: (1) local buckling analysis, and (2) global buckling analysis (general panel instability). The follow- 
ing sections describe these approaches. 


Local Buckling 


The purpose of local buckling analysis was to study the buckling behavior of a local weak region of 
the panel. This weak region is identified as a rectangular flat plate region bounded by two legs of the re- 
inforcing hat located at the center of the global panel (left diagram of fig. 2). The analysis looked at the 
buckling behavior of this rectangular flat plate (slender strip). Because the reinforcing hat has high flexural 
rigidity, the four edges of the rectangular plate were assumed to be simply supported (right diagram of 
fig. 2). From reference 2, the shear buckling stress (x xy ) cr in the rectangular flat plate may be written as 
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Figure 2. Shear buckling of hat-stiffened panel analyzed using a simple model. 


where k is the shear buckling load factor, which is a function of panel aspect ratio — . For — = 1 (a 
y b a a 

square panel), k = 9.34 ; and for — = 0 (an infinitely long panel), Jt = 5.35 . For intermediate val- 
b y a } 

ues of — , k may be found from a parabolic curve-fitting equation of the form (ref. 2) 
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The curve described by equation (2) is shown on the left in figure 3. The panel shear load TV for the hat- 
stiffened strip (left diagram of fig. 2) may be written in terms of shear flows (fig. 4) as (ref. 4) 
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( 3 ) 


N xy = 4 , +q c 


where q x and q c are, respectively, the shear flows in the flat panel and the hat, and are given by (ref. 4) 
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where 4^- is the panel twist. 


dxdy 

From equations (4) and (5), the ratio q ( iq c may be calculated. Then, from equation (3), the panel 
shear buckling load ( N ) cr of the hat-stiffened strip may be calculated as a function of ( T xy ) ^ (eq. ( 1 )). 


Global Buckling 


In the global buckling analysis (general instability analysis), the complex panel was represented by a 
homogeneous anisotropic panel having effective elastic constants. These effective elastic constants must 
be calculated first (ref. 3). This analysis is similar to the conventional buckling analysis of a sandwich 
panel with one face sheet removed. 

By using the small-deflection theory developed for flat sandwich plates (ref. 5) and solving the shear 
buckling problem of the hat-stiffened plate using the Rayleigh-Ritz method of minimizing the total poten- 
tial energy of a structural system (refs. 5 through 9), the following shear buckling equation is obtained: 
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Figure 5. Segment of hat-stiffened flat panel. 


where 
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and / c is the moment of inertia, per unit width, taken with respect to the panel neutral axis T| (fig. 5), 
and is given by 
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where the nondimensional coefficient S is defined as 
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where the nondimensional parameters D , D , and D are defined as 


£> f = 


“ 3\F 


/_ r 


2'^ 3 cos 2e + ?-£ 
3/ / 


w-(*n 


/? 

+ *- 


2 ( ^) 2 9 _ ( 1 - cos 0 ) + ^ (0- sin0cos0) 


.Rb 

h 


R \ 2 


h 2 t L 

C C 


2^sin 2 0 + £ (0 - sin 0 cos 0) 

K h c J 


(29) 


(30) 


+ ^ |(^-)0-2^|(0- sin0) -jj- ( 1 -cos 0) j^l- ^ ( 1 -cos 0) j 

—( 2 £ sin0cos0 + £sin 2 0^ 

i . 2 * \ K K ) 


(31) 


where 


D h = \ ( 3 sin 2 0 + \ ( £-0 + \{-y\ 

y 3 \hj 2 \h c 2 h c I f l 

- ^2-3^ (0-sin0) + ^sin0 (1 - cos0) 




■f-— + 2j-cos 2 0 + £ (0 + sin0 cos 0) 


L h c<f 


h 


f = |t/,+/ 2 ), * = | [p-|(/, + /2>]. >c - n'"’ '/ = ^ 


1 , 3 

5 


(32) 


(33) 


The shear buckling equation (6) yields a set of homogeneous equations associated with different values of 
m and n. This set of equations may be divided into two groups that are independent of each other: one 
group in which m ± n is odd (that is, antisymmetrical buckling), and the other in which m ± n is even (that 
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For m ±n = odd (antisymmetrical buckling): 
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where the nonzero off-diagonal terms satisfy the conditions m ^ i,n kj, m ± i= odd, and n ± j = odd. 

Notice that the diagonal terms in equations (34) and (35) came from the first term of equation (6), and 
the series term of equation (6) gives the off-diagonal terms of the matrices. The 12 x 12 determinant was 
found to give sufficiently accurate eigenvalue solutions. 
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Figure 6. Buckling shape of hat-stiffened panel under shear loading (finite element analysis by 
W. Percy, Mc-Donnell-Douglas; full-panel model). 

This local shear buckling load prediction is slightly higher than the value (N xv ) cr = 900 lb/in calcu- 
lated from finite element buckling analysis carried out by W. Percy of McDonnell-Douglas. Figure 6 
shows the shear buckling shape of the hat-stiffened panel based on Percy’s full-panel finite element model. 
Clearly, the panel is under local buckling rather than general instability. The local buckling analysis pre- 
dicts a slightly higher value of (N ) cr because the four edges of the rectangular plate strip analyzed were 
assumed to be simply supported. In reality, those four edges are elastically supported. 

Global Buckling 

To find the order of the determinant (review eqs. (34) and (35)) for converged eigenvalue solutions, 
several different orders of the determinants were used for the calculations of k . The eigenvalues were 
found to have sufficiently converged beyond order 10. In the actual calculations of k x , the orders of the 
determinants were taken to be 12, which were shown in equations (34) and (35). The eigenvalue solutions 
thus obtained give the following lowest values of k : 

xy 

m±n = even: k =1.89 (41) 

m±n = odd: k = 1.93 (42) 

Thus, the square panel will buckle symmetrically. Using k = 1.89, the panel shear buckling load ( N ) 
may be calculated from equation (10) as ' cr 

(N xy ) cr = 4, 296 lb/in (43) 

Personal communication with author. 
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